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ABSTRACT 

NASA is developing a 15-kW Hall thruster to support future NASA missions. This activity is 
funded under the Space Technology Mission Directorate Solar Electric Propulsion Technology 
Demonstration project. As a part of the development process, the far-field plume characteristics of the 
NASA-300M, a 20-kW Hall thruster, were studied. The results will be used to study how various aspects 
of the operation of this thruster affect the overall performance. This data will be used to guide future 
design work and serve as a baseline for comparison to a magnetically shielded version of the NASA- 
300M that will be tested in the future. For this study, a far-field Faraday probe was swept in a polar 
fashion to map the ion current density. An ExB probe (Wien filter), two retarding potential analyzers, and 
a Langmuir probe were mounted at a fixed location on the thruster axis in the far-field plume. The data 
reduction method followed recommendations in recent studies by Brown, Reid, and Shastry with 
modifications that are tailored to the plasma plume environment of high-power Hall thrusters. Results from 
this and prior testing show that the plume is richer in doubly-charged ions, larger in spatial extent, and 
capable of greater probe heating than lower power thrusters. These characteristics require special 
treatment in experimental setup and data analysis, which will be the main focus of this paper. In 
particular, covered topics will include a new, more accurate, method of integration for analysis of the ExB 
probe data and effect of secondary electron emission on the Faraday probe data. 
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INTRODUCTION 



In 2010, NASA established the Human Exploration Framework Team (HEFT) to analyze 
exploration and technology concepts and provide inputs to the agency's senior leadership on the key 
components of a safe, sustainable, affordable, and credible future human space exploration endeavor for 
the nation. 1 The team concluded, in part, that the use of a high power (i.e. on the order of 300 kW) solar 
electric propulsion system could significantly reduce the number of heavy lift launch vehicles required for 
a human mission to a near earth asteroid. 1,2 Hall thrusters were found to be ideal for such applications 
because of their high power processing capabilities and their efficient operation at moderate specific 
impulses, which leads to reduced trip times that are necessary for such missions. 2 ' 4 Recent electric 
propulsion system model estimates that considered factors such as cost, mass, fault-tolerance, cost 
uncertainty, complexity, ground test vacuum facility limitations, previously demonstrated power 
capabilities, and possible technology limitations have shown that Hall thrusters operating at power levels 
of 20-50 kW are strong candidates for human exploration missions operating at total powers up to 500 
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As an intermediate step towards the aforementioned higher power Hall thruster, the Space 
Technology Mission Directorate (STMD) is investing in the development of a 15-kW Hall thruster. This 
thruster will serve as a stepping stone to higher power and provide a state-of-the-art propulsion system 
for the Solar Electric Propulsion (SEP) Technology Demonstration Mission (TDM). SEP TDM, initially 
announced in 2011, is aimed at demonstrating new cutting edge technology in solar arrays, spacecraft 
bus, and electric propulsion that will increase the maturity of these SEP technologies for future 
commercial use. Once these technologies have been demonstrated, they are expected to enable higher 
performance low-Earth orbit to geosynchronous Earth orbit (LEO to GEO) transfers as well as a number 
of other near-Earth orbit transfers and station-keeping maneuvers. In longer terms, these technologies 
will reduce mission costs for NASA interplanetary missions and will serve as a precursor to higher power 
systems for manned interplanetary exploration. 

The NASA Glenn Research Center (GRC) is partnering with the Jet Propulsion Laboratory (JPL) 
to study high power Hall thrusters and obtain data for designing the aforementioned 15-kW Hall thruster. 
One of the activities prior to the design and fabrication of the 15-kW Hall thruster will be a set of studies 
comparing the NASA-300M, a 20-kW laboratory Hall thruster, with a magnetically-shielded version of the 
same thruster called the NASA-300MS. The NASA-300M was developed from 2004 to 2005 for project 
Prometheus but was not tested until recently due to project cancellation. 6 With the recent renewal in 
interest towards high-power Hall thrusters, the performance of the NASA-300M was characterized in 
201 0. 6 The study that will be described in this paper was a far-field probe interrogation of the 300M 
plasma plume. Information obtained from this study will form the baseline for comparison with the 300MS 
in a future study. These studies combined with the modeling effort by JPL will help confirm the principles 
of magnetic shielding when scaled up to power levels above 10 kW. 

For this study, a far-field Faraday probe was swept in a polar fashion to map the ion current 
density. An electric-field-cross-magnetic-field (ExB) probe, also known as a Wien filter, two retarding 
potential analyzers (RPAs), and a Langmuir probe were mounted at a fixed location on the thruster axis in 
the far-field plume. The ExB probe was used to measure species fraction of the multiply-charged xenon 
species as compared to singly-charged xenon; the RPAs were used to measure energy per charge 
information of the plume; and the Langmuir probe was used to correct the RPA data. One of the RPA is 
newly designed by the Air Force Research Laboratory (AFRL) and was being validated. 


The main focus of this paper will be on the experimental setup and modifications to traditional 
data analysis methods that came about as a result of special characteristics of high-power Hall thruster 
plumes. The data reduction methods are based on recent studies by Brown, 7 Reid, and Shastry 9 with 
modifications that are tailored to the plasma plume environment of high-power Hall thrusters. Results from 
this and prior testing show that the plume is richer in doubly-charged ions, larger in spatial extent, and 
capable of greater probe heating than lower power thrusters. These characteristics require special 
treatment in experimental setup and data analysis, which will be the main focus of this paper. In 
particular, covered topics will include a new, more accurate, method of integration for analysis of the ExB 
probe data and effect of secondary electron emission on the Faraday probe data. The test methodology 
and data analysis carried out for this paper draws heavily from the knowledge gained during those tests. 


II. EXPERIMENTAL SETUP 


For the sake of brevity and to simplify plot labeling, NASA-300M will sometimes be referred to as 
the 300M in this paper. Operating conditions are labeled as www V, yy.y A, where www is the discharge 
voltage in volts and yy.y is the discharge current in amperes. Unless otherwise noted, all spatial positions 
presented in this paper have been normalized by the mean thruster diameter of the thruster. Mean 
thruster diameter (MTD) is defined as the average of the inner and outer discharge-channel wall 
diameters. 
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Figure 1. NASA-300M mounted on a thrust 
stand. 

A. Thruster, Thrust Stand, and Test Matrix 



Discharge current, A 

Figure 2. Operating conditions tested for this 
study. 


The NASA-300M is a magnetic-layer Hall thruster. It was designed based on the principles 
outlined in Manzella’s dissertation work. 10 A photograph of the 300M mounted on the thrust stand is 
shown in Fig. 1. This thruster was originally developed for high-specific-impulse missions 6 and has lens- 
type magnetic field topology. 11,12 The thruster ran with a centrally-mounted hollow cathode derived from 
the discharge cathode for NASA’s Evolutionary Xenon Thruster (NEXT). 6 For the data presented in this 
paper, the magnetic field settings of the thruster were optimized to give a shape roughly symmetric about 
the channel centerline while maximizing anode efficiency. The cathode flow fraction was 8%. 

The nominal operating condition of NASA-300M is 500 V, 40 A. For this test, the 300M is 
operated at 10 conditions spanning discharge voltages of 200 to 500 V, discharge currents of 20 to 67 A, 
and discharge powers of 10 to 20 kW. The anode efficiency of this thruster, defined in Eq. (1), varies from 
57 to 69% over the tested operating conditions. 


In this equation, T is thrust, rh a is anode mass flow rate, l d is discharge current, and V d is discharge 
voltage. The boron nitride channel walls of the 300M have accumulated <110 hours of operation time 
prior to the start of this study. 

The thrust stand used in this study is an inverted pendulum thrust stand designed by Haag. 13 It is 
actively cooled during operation. The nominal accuracy of this thrust stand is ±2%. Long-term thermal 
drift is corrected by measuring the thrust signal with all gas flow to the thruster off and then assuming a 
linear change in the zero-thrust value over time. The maximum thermal drift was found to be no greater 
than 9 mN. 

The throttling table for this study was designed to target a variety of discharge voltages in order to 
characterize the thruster for both high thrust-to-power (low specific impulse) and high specific impulse 
operation. Many of the conditions chosen also enable constant-current and constant-power comparisons. 
Figure 2 shows a diagram of all operating conditions tested for the present study of the 300M. 

B. Vacuum Facility and Motion Stages 


Testing was conducted in NASA GRC’s Vacuum Facility 5. This cylindrical facility is 4.6 m in 
diameter and 18.3 m long and is pumped with a set of cryo-panels and 20 oil diffusion pumps. The 300M 
was mounted on a thrust stand in the main volume of the vacuum facility to ensure low pressures near 
the thruster during operation. Facility pressures were monitored with one Pirani gauge and six ion 
gauges, one of which is mounted next to the thrust stand. Facility pressure uncertainties were estimated 
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Figure 3. Photograph of the diagnostics setup. 

by the manufacturer to be within ±6% of reading. The pressure reading, corrected for xenon, next to the 
thrust stand was less than 4.0x1 0" 3 Pa (3.0x1 O' 5 Torr) throughout testing. 

A gas feed system made of commercially available components was used to feed xenon to the 
thruster. A 1000-sccm mass flow controller provides propellant to the anode while a 200-sccm controller 
provides propellant to the cathode. Both controllers were calibrated using xenon prior to testing. Typical 
uncertainties reported herein were ±1.0% of reading for the anode mass flow and ±1.1% of reading for 
total mass flows. 

The commercial discharge power supply used is capable of supplying 2000 V at over 100 A with 
a 5.5 mF output filter capacitance. A 15.3 mF capacitor bank was connected in parallel with the power 
supply at the electrical feed-throughs of the vacuum facility to filter discharge current oscillations. In 
addition, separate commercial power supplies were used to power the cathode heater, cathode keeper, 
and electromagnets. 

A farfield Faraday probe was mounted onto a commercially available three-axis belt-driven 
motion system. The motion system provides 2D rectilinear motion and probe rotation. Positioning 
accuracy is ~1 mm and <0.1°, respectively. Due to time constraints, data were only taken with the 
Faraday probe at 5 MTD away from the center of the thruster at exit plane. An ExB probe, two RPAs, and 
a Langmuir probe were mounted to a probe tower with accompanying shielding and shutters protecting 
the ExB probe and the RPAs. This farfield probe tower was attached to a vertical motion stage located 
~27 MTD downstream of the thruster exit plane and positioned behind a Grafoil shield. During the data 
acquisition process, each probe was moved vertically to a point roughly along the central axis of the 
thruster before measurements were taken. The probe tower retreated behind the shield once data 
acquisition was completed. Figure 3 shows a photograph of the diagnostics setup. Figure 4 shows a 
photograph of the farfield probe tower. 

C. Description of Probes and Data Acquisition System 

This section will describe the probes used in this study in the following order, Langmuir probe, 
RPA, ExB probe, and Faraday probe. Unless otherwise stated, all probe biases were applied with 
commercial power supplies. The RPAs and ExB probe had leaf-blade shutters in front of them. 

The Langmuir probe consists of a single tungsten wire protruding from an aluminum oxide tube. 
This probe was used to obtain the plasma potential at the same location as the RPA so that the RPA data 
can be corrected for this potential. This Langmuir probe was swept at 10 Hz for 1 second at each test 



Figure 4. Photograph of the 
farfield probe tower. 
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point. The probe was connected to a custom circuit box where the probe current was passed through a 
shunt and the signal fed to an isolation amplifier. The bias voltage was passed through a voltage divider 
and fed to another isolation amplifier. Signal from the amplifiers were fed to a NI-9205 data acquisition 
module attached to a Nl cDAQ-9178 unit. 

Two four-grid RPAs were used in this study. These RPAs were used to approximate the average 
ion energy per charge. RPA 1 is of GRC design and was found based on prior experience to not work 
well at closer than -25 MTD for medium to high power Hall thruster testing. For this reason, the probe 
tower was located at -27 MTD from the thruster. RPA 2 is of AFRL design, courtesy of Daniel L. Brown 
and Joseph Blakely. For each RPA, the electron suppression and repelling grids were biased to -30 V 
with respect to facility ground while the ion retarding grid voltage was swept. The ion retarding grid was 
biased by a Keithley 2410 sourcemeter while the collected current was measured by a Keithley 6485 
picoammeter. Results collected during the study show the two RPAs agree with each other and RPA 2 
shows superior noise characteristics and stand-off against break down. Presented results will all be from 
RPA 2. 


The ExB probe is a commercial product of Plasma Controls, LLC. It was used to measure 
charged species current fraction. The ExB probe weighs -2 kg and is -270 mm long. This ExB probe was 
the result of a Small Business Innovation Research project and has a proven history of usage. 14 Ions 
entering this ExB probe must first travel through a collimator, which ensures that the ions enter the ExB 
filter section along the probe axis. Perpendicular electric and magnetic fields then form the filter section 
where only charged particles of the right velocity can pass through to the electron suppression plate. A 
drift section separates the filter section and the electron suppression plate. The drift section enhances 
probe resolution and makes sure bias on the suppression plate does not affect the filter section. The 
electron suppression plate has a small orifice through which ions can pass through. It was biased at -30 V 
with respect to facility ground to suppress secondary electron emission from the collector. Ions passing 
through the suppression plate were collected by the collector. The collimator used in this study has an 
entrance and exit orifice diameter of 0.76 mm. The bias plates that form the electric field in the filter 
section have a gap distance of 3 mm. The electron suppression plate orifice diameter is -1.5 mm. The 
collector is large enough to collect all ions that pass through the electron suppression plate orifice. A 
thermocouple was attached to the case of this ExB probe to make sure it did not overheat and de- 
magnetize. During testing, the main bias plate voltage was sw< 
picoammeter, which also measured the collector current. 

The Faraday probe, which is of GRC design, was used to 
measure ion current density. Figure 5 shows a cross-sectional diagram 
of this Faraday probe. The collector area is a circle of diameter 17.39 ± 

0.02 mm. The front area of the collector is 238 mm 2 . The probe has an 
overall diameter of -28 mm. The collector and guard ring are made of 
molybdenum (to minimize secondary electron emission) and the 
insulating back is made of Macor. The collector and guard ring 
connections are mated to wires behind the insulating back and covered 
with fiberglass tape. The Faraday probe was set at 5 MTD away from 
the thruster because it needed to be far enough away from the thruster for the source of the thruster 
exhaust to be considered a point source. At the same time, the probe cannot be set too far away due to 
the effects of charge-exchange (CEX) between plume ions and background neutrals. During testing, the 
probe was set to different bias voltages (-20 to -40 V), and -30 V with respect to facility ground was 
selected to ensure electrons were properly repelled but sheath expansion was negligible. Collected 
current was measured via a shunt and an isolation amplifier. The output of the amplifier was fed to 
another NI-9205 data acquisition module attached to the aforementioned Nl cDAQ-9178 unit. 

The experiment was conducted via a dedicated data acquisition computer, which was running a 
custom Labview program. The motion stages were controlled by a motion controller that received 
movement commands from the computer. Encoder signals from the stages were fed into an encoder 
controller, which reports to the computer. The Nl cDAQ-9178 data acquisition device and the various 
Keithley devices fed data to and were commanded by the computer. During the experiment, the computer 
automatically orchestrated the sequence of events that activated the various motion stages, shutters, and 
probes. 


;pt biased by a Keithley 6487 
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Figure 5. Diagram of the 
farfield Faraday probe. 
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III. DATA REDUCTION 


Data reduction was carried out in a specific sequence because, for the most part, each step of 
the analysis depended on the data obtained in the preceding steps. Though some of the analysis steps 
can be carried out independently, the accuracy of the result is decreased. For this study, the data 
analysis sequence is as follows, Langmuir probe analysis, RPA analysis, ExB probe analysis, and 
Faraday probe analysis. Additionally, pressure data is analyzed to help correct the ExB probe data for 
facility-background CEX effects. Differences between the presented data reduction procedures and 
traditional approaches will be highlighted with explanations for the rationale behind the new procedures. 

A. Langmuir Probe Analysis 

Langmuir probe analysis was carried out using simple cylindrical Langmuir probe theory. 15 Since 
the Langmuir probe was swept at 10 Hz for 1 second, 20 ramps were present in each data file. These 
data were averaged then smoothed using a Savitzky-Golay filter. 16,17 Then, the derivative of the probe 
current with respect to the probe voltage was plotted against the probe voltage. The voltage 
corresponding to the highest amplitude value on this plot is roughly equal to the plasma potential. This 
potential is needed to correct the RPA reading because the RPA ion retarding grid was biased with 
respect to the facility ground. The true filter voltage was equal to the ion retarding grid bias voltage minus 
the local plasma potential. 

B. RPA Analysis 

RPA analysis was carried out by first smoothing the RPA trace then taking the negative of the 
derivative of the collector current with respect to the bias voltage on the ion retarding grid. The result, 
plotted against the bias voltage, is proportional to the ion energy per charge distribution function. 18 
Traditionally, the most probable voltage is taken as the average ion energy per charge. The most 
probable voltage is defined as the voltage where the amplitude of the plot maximizes. However, due to a 
combination of noise and broad ion distribution, the data taken for the 200 V, 49.8 A case had a small 
peak on top of the main population that skewed the location of the most probable voltage. To deal with 
this issue, the average ion energy per charge was calculated via ensemble-averaging using only the part 
of the trace where the amplitude exceeded half of the maximum amplitude. This averaging approach will 
be referred to as the threshold-based averaging approach with a 50% threshold. For clean traces, this 
approach yield values that are at most 3-4 V different from the most probable voltage. For noisy traces, 
this approach yield values that are insensitive to noise spikes. Figure 6 shows plots from applying the 
threshold-based averaging approach to the RPA data obtained for the 200 V, 49.8 A condition. The black 
dashed vertical line indicates the location of the most 
probable voltage, the red solid vertical line indicates 
the result of using the threshold-based averaging < 
approach with the 50% threshold, and the red dashed = 
horizontal line indicates the 50% of maximum | 
threshold. < 

CU 

Note that in theory the most accurate result is ~ 
obtained by ensemble-averaging the entire RPA trace. 

However, doing so often produces unphysical results 
because the ion energy per charge distribution as 
measured by the RPA is typically much broader than = 
the real distribution due to the wide acceptance angle 
of the RPA. This is also the reason the traditional g 

approach relies on the most probable voltage. Using „ 

the 50% threshold-based ensemble averaging 1 
essentially strikes a balance between excluding JP 
broadened data and maintaining some degree of 
noise insensitivity. 
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6. Sample RPA analysis plots. 
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The ion energy per charge obtained from the RPA, corrected by the Langmuir probe 
measurement, was inputted into the CEX model used in the analysis of the ExB probe data. CEX 
depletes lower-charge-state species more than higher-charge-state species, so the effect of ion CEX 
collisions with facility background neutrals must be accounted for. This effect is especially prominent for 
high-power thruster tests. One can approximate the ion energy per charge as the discharge voltage. 
Doing so removes the dependence of the ExB probe analysis on the RPA and Langmuir probe data, 
though with a small increase in uncertainty. See Shastry’s work for an estimate of this uncertainty. 9 


C. Traditional Approaches for ExB Probe Analysis 


The ExB probe is strictly a velocity filter. Since different charged species are accelerated to 
different velocities by going through roughly the same potential drop, they will show up as different peaks 
when interrogated by the ExB probe. Traditionally, the simplest ExB probe analysis approach is to take 
the height of each peak and assume they are proportional to the current of the respective species. This 
approach is identical to assuming mono-energetic particles for each species. 19,20 The peak height 
approach is also used as a last resort if the width of the ion velocity distribution measured by the ExB 
probe is around the same or less than the ExB probe’s velocity resolution. In such a case, the large ExB 
probe velocity resolution causes artificial broadening of the recorded traces. 

If the ExB probe velocity resolution is at least several times smaller than the width of the 
measured ion distribution, the preferred method for analyzing ExB probe data is by integration. Most 
integration schemes involve fitting some general distribution equation form to the ExB probe trace. Then, 
the area under the curve is calculated either by geometry or by numerical integration. Curve-fit and 
integration are used together because the ion distributions from different charged species tend to overlap 
in the regions between peaks, making integration of only one species at a time impossible without some 
kind of assumption on the form of the distribution. The domain of the curve-fit for each species generally 
need to be bounded such that only current signal generated by one species is being fitted at a time. 

One of the earliest examples of the traditional integration approaches is the half-width-at-half- 
maximum (HWHM) triangle fit, proposed by Beal. 21 In this approach, the current due to each species is 
assumed to be proportional to the product of the peak height and the HWHM value. This is equivalent to 
fitting a triangle to the measured ExB probe trace for each species in such a way that the peak and 
HWHM value of the triangle match those of the species’ distribution. 

Another curve-fit approach is the Gaussian fit, proposed by Linnell. 22 In this approach, each peak 
in the ExB probe trace is fitted to a Gaussian curve within a limited curve-fit domain. The curve-fit domain 
is individually selected for each species so that most of the current signal inside the domain is guaranteed 
to come from only the associated species. The integrated area under the fitted curve is assumed to be 
proportional to the current fraction of the associated species. 

A third curve-fit form is the variable exponent fit proposed by Kim. 23 Kim attempted to find an 
exponent-based equation form that best fit the ExB probe data of the SPT-100 and found that the 
distribution function typically lies between the Gaussian, which scales with exp(-u 2 ), and the Druyvesteyn, 
which scales with exp(-u 4 ), where u is the particle velocity. Since the Gaussian function describes a 
homogeneous fluid at thermal equilibrium and a Druyvesteyn function describes a charged population in a 
uniform electric field with elastic collisions against a uniform background, the rationale is that the ion 
distribution of the Hall thruster exhaust must lie somewhere in between. Thus, Kim derived a formula 
where the exponent is left as a fitting parameter. Equations for the various fitting equations will be shown 
in the next section. Note that Kim’s work is the first published instance where the ExB probe was used to 
study a Hall thruster and strongly influenced all subsequent work in this area including all other cited work 
on ExB probe studies. 

Once curve-fitting is done, current fractions can be calculated via Eq. (2), 9 
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where Qj is the current fraction of the j-th species, lj is the current of the j-th species, n f is the density of 
the j-th species, and Z t is the charge of the j-th species. For the traditional integration approaches, 
regardless of the form of the fitting function, the current of the j-th species is calculated from the total area 
under the curve found in the curve-fitting step multiplied by constants that are the same across all 
species. Since these constants are the same, the current fraction can be calculated simply as the ratio of 
the integrated probe current for a given species divided by the sum of integrated probe current. The 
species fraction can be calculated using Eq. (3), 9 

_nj_ fyz j 3 ' 2 

1 ' Z" k - S^/Zk 3 ' 2 < 3 > 

k k 

where Q is the species fraction of the j-th species. Note, that the last term of Eqs. (2) and (3) only applies 
if one can assume negligible variation in acceleration voltage for different species. The effect of 
background CEX is typically corrected at this step by modifying the integrated current and integrated 
density. The CEX correction scheme used in this study is described in Shastry’s work 9 and reproduced in 
Eqs. (4)-(6), 


(J/JoW = ex P(- n 0°1 z )> 

o-| =87.3-13.6 logM) 

(4) 

(J/J 0 ) Xe 2 + = exp(-n 0 a 2 z), 

a 2 =45.7-8.9 log(2V 2 ) 

(5) 

(J/J 0 ) Xe 3+ = exp(-n 0 a 3 z), 

o 3 =16.9-3.0 log(3V 3 ) 

(6) 


where J is the current density at distance z away from the thruster exit plane, J 0 is the current density at 
the thruster exit plane, and n 0 is the average background neutral density. The numerical formulas for the 
CEX cross-sections, a 1; a 2 , and o 3 , are given in units of A 2 (10" 2 ° m 2 ). The ion energy per charge, V l5 V 2 , 
and V 3 for Xe + , Xe 2+ , and Xe 3+ , respectively, are assumed to be equal to the average ion energy per 
charge as measured by the RPA. The value of n 0 is calculated by converting the average of the pressure 
measurement from four ion gauges. Two of these gauges were located on the chamber wall at the same 
axial position as the farfield probe tower, one of the gauges was next to the thrust stand, and the last of 
the gauges was on the chamber wall at an axial location close to the thruster. 

In summary, the typical approach is to perform a curve-fit for each peak using the VDF form of 
choice, integrate the current, apply the CEX correction, calculate the current fraction, and, if necessary, 
calculate the species fraction from the current fraction. 9 The current fraction can then be used in any of a 
number of phenomenological efficiency models. These models are designed to help researchers identify 
trends in the various physical processes that drive the overall efficiency of a Hall thruster. The 
phenomenological efficiency model used in this study will be described in a later section. 

D. New Approach for ExB Probe Analysis 

The approach used to analyze the ExB probe data for this study was derived from scratch and 
closely matches the derivation by Kim. 24 However, Kim did not publish the exact integration formula used 
to calculate the species fraction and did not try to calculate the current fraction (that was not his 
objective). To the best of the authors’ knowledge, the integration formulas presented in this paper are 
original and not published in any prior work. 

To figure out the relationship between the current fraction and the ExB probe trace we will first 
look at the physical description of what exactly the ExB probe is collecting at a given plate bias voltage. 
This is shown in Eqs. (7), (8), and (9), 

lp(V P ) = ^l Pjj (Vp) 
j 

fUn+AU 

l Pi (Vp) = Z i eApn i F uf,(u)du 

’ J J Ju 0 -Au ‘ 


( 7 ) 

( 8 ) 
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( 9 ) 


Ep V P 

B 0 B 0 D f 


where l P is ExB probe collector current, V P is the plate bias voltage, l Pij is the contribution to the measured 
current by the j-th species, e is the elementary charge, A P is the probe entrance area, u is particle 
velocity, u 0 is the particle velocity corresponding to the applied electric and magnetic field magnitude, E 0 
and B 0 , respectively, as shown in Eq. (9), AUj is the half-width of the ExB velocity resolution (its charge 
dependence is derived in the appendix), fj is the velocity distribution function (VDF) of the j-th species, 
and D f is the gap between the ExB bias plates. 

For an ExB probe with good velocity resolution, that is to say AUj « u 0 , the integral term can be 
approximated as fj(u 0 )*2AUj. The derivation of the dependence of AUj on other basic parameters is shown 
in the appendix. The final result is replicated in Eq. (10) for convenience, 


A Uj =^G 
J Z;eB n 


( 10 ) 


where m is the mass of the particle and G is a function of only geometric constants. Taking the AUj « u 0 
approximation and substituting Eq. (10) into Eq. (8), we arrive at Eq. (11), 


l PJ (Vp) = K 0 Uo 3 n j f j (uo), 


( 11 ) 


where K 0 is a constant for a given probe design and propellant atomic mass. Further substituting Eq. (9) 
into Eq. (11), we arrive at Eq. (12), 

lp, j (Vp) = K 1 n j Vp 3 f j (K 2 V P ), K 1S — , K 2 = — L- (12) 

J J J (B 0 D f ) 3 B 0 D f v ' 


where Ki and K 2 are constants for a given probe design. Depending on which VDF form is in use, the 
associated equation is substituted into fj, and Eq. (12) is fitted to the ExB probe trace to solve for any 
fitting constant. The formulas presented so far, though appearing in different forms, matches previously 
published work. 9,23 The formulas to follow are original and diverge from past work. 

To calculate the species fraction and current fraction, we now need the integral equations relating 
each quantity to the ExB probe current. To find the species fraction, we begin by re-arranging the V P 
term, then integrating both sides of Eq. (12) with respect to V P from 0 to °°. Since n f does not depend on 
V P it can go outside the integral. The result is shown in Eq. (13). Since the integral of fj from 0 to 00 is 
unity, we can further simplify the result, which is shown in Eq. (14). 


K,n,jl i (K 2 Vp)dVp=J^^dV l 


n i = 

J K 


1 r c 

k7 Jo 


'lpj(Vp) 


dVp 


V P 


(13) 

(14) 


Equation (14) relates nj to the ExB probe current and can be used to calculate the species fraction. The 
integration formula for species fraction is shown in Eq. (15). 



k 


r 00 lp,j(Vp) 

t) v P 3 


dV P 


y f°° !p,k( v p) 
k Jo Vp 3 


dVp 


(15) 
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To find the current fraction, we begin by relating the current of the j-th species, lj, to the ExB probe 
current. The current of the j-th species is defined in Eq. (16). Substituting Eq. (9) into (16) yields Eq. (17). 


poo 

lj =Z j eApn j J o U 0 fj(Uo)dUo 

(16) 

= ZjeApK 2 2 nj f°V P fj(K 2 Vp)dV P 
JO 

(17) 


Next, we rearrange Eq. (12) by dividing both sides by V P , then integrate both sides with respect to V P from 
0 to 00 to obtain Eq. (18). Substituting Eq. (18) into Eq. (17) yields Eq. (19). 


nj f°Vp fj(K 2 Vp)dV P ^r^dV, 
Jo Kt Jo y P 2 


ZjeA P K 2 <*colpj(Vp) 
Ki Jo v P 2 


-dV P 


(18) 

(19) 


Equation (19) relates lj to the ExB probe current and can be used to calculate the current fraction. The 
integration formula for the current fraction is shown in Eq. (20). 


Qj = 


5 > 



lp,j(Vp) 

Vp 2 


dVp 



lp,k(Vp) 

Vp 2 


dVp 


( 20 ) 


In the new ExB probe analysis approach, the process starts with selecting a form for the VDF, then 
inserting the form into Eq. (12) and performing a curve-fit to the ExB probe current as a function of bias 
plate voltage for each species. Using the result, Eq. (15) and (20) are then numerically integrated to 
obtain the species fraction and current fraction, respectively. The CEX correction is applied prior to 
normalization. For all ExB probe analysis carried out in this study, the curve-fit process starts by curve- 
fitting the highest charge peak present (Xe 4+ ), subtracting it out from the probe trace, then proceeding to 
lower charge peaks. Xe 4+ data is not used in the final current and species fraction calculation because it 
typically comprises less than 0.5% of the beam current. Flowever, Xe 4+ curve-fit is still performed and 
result subtracted to make curve-fitting of lower charge states easier. 

Listed in Eqs. (21)-(24) below are the forms of several VDFs that were tried in the current study. 
They are, in order, triangle, Gaussian, variable exponent, and twin Gaussian. 


Triangle: 

Gaussian: 
Variable Exponent: 
Twin Gaussian: 


, . . _ J(-a/b) * |u — c| + a, c-b<u<c + b 
j U { 0, else 

fj(u) = a*exp[-b(u-c) 2 ] 
fj(u) = a*exp[-b|u-c| n ] 

fj(u) = a-| *exp[-b-| (u-c-|) 2 ] + a 2 *exp[-b 2 (u-c 2 ) 2 ] 


( 21 ) 

( 22 ) 

(23) 

(24) 


where a, b, c, n, a-i, bi, Ci, a 2 , b 2 , c 2 , are fitting constants. Note the triangle form listed here is not quite 
equivalent to the FIWFIM triangle method described earlier because this VDF will actually be fitted to the 
ExB probe data as opposed to having only its peak and width matched. The triangle form is tried because 
it is a close match to the FIWFIM triangle method and because the FIWHM triangle method is incompatible 
with the new integration formulas. The twin Gaussian form was tried in an attempt to capture the tail that 
is often observed in the VDFs of accelerated ions in a Flail thruster. While this form generally gives a 
much better fit to the data than any other form, it can generate an erroneous result if the curve-fit 
boundary is not properly chosen. For instance, if the boundary includes peaks from two different species, 
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the computer will fit the twin Gaussian form to both peaks and treat them as a single species. The 
operator must study each curve-fit plot to ensure this is not occurring during data analysis. 

As a side note, the traditional integration formula for the current fraction can be written as Eq. 

(25). 


/•oo 

I, J 0 'p,|(Vp)dV P 
T [ ip t: (Vp)dVp 

k 

k 


(25) 


While Eq. (25) is mathematically incorrect, it turns out to be a pretty good approximation of Eq. (20) for 
analyzing the plume of a Hall thruster, especially when the VDF is narrow. A comparison between the 
traditional and the new integration formulas will be made in the results section. 


E. Faraday Probe Analysis 


Faraday probe data are used to calculate the plume divergence angle and the total ion beam 
current. The cosine of the momentum-weighted plume divergence angle is defined as the average axial 
velocity of the particles divided by the average velocity of the particles. However, momentum-weighted 
divergence angle is difficult to measure. The typical approach is to measure the charge-weighted 
divergence angle, which is approximately equal to the momentum-weighted divergence angle if the 
current fraction is roughly constant across the interrogated domain. For a polarly-swept probe, Eq. (26) 
can be used to calculate the charge-weighted divergence angle, 

9 p7l/2 

27rRpp I j(0) cos0sin6d0 

< cos 5 >= 0 W2 (26) 

2n Rpp 2 J j(0) sin0d0 

where 5 is the charged-weighted divergence angle, 0 is the polar angle and is equal to 0° for particles 
traveling parallel to the firing axis, and j(0) is the ion current density as a function of the polar angle. R F p is 
the distance from the Faraday probe collector to the thruster center at exit plane and is constant for a 
polar sweep. Note the denominator is equal to the total ion beam current. 

For a nude Faraday probe with a guard ring like the one used in this study, the effective collector 
area is not exactly equal to the collector frontal surface area. Current that enters the gap between the 
collector and the guard ring can be collected by the side surfaces of the collector. 25 According to work by 
Brown, the current entering the gap is collected by the collector and the guard ring in ratio proportional to 
the ratio of exposed gap area. 25 For the probe design used in the current study, the area inside the gap is 
dominated by area connected to the guard ring (See Fig. 5). However, there is enough area connected to 
the collector that some level of correction is needed. Note that only the part of the gap with direct 
exposure to the incoming ion beam is used in the gap area calculation. Using the approach 
recommended by Brown, the effective collection area is ~4% greater than the collector frontal area. The 
effective collection area is used for all Faraday probe analysis. 

Much like in the analysis of the ExB probe data, CEX is a factor in Faraday probe data analysis. 
However, whereas CEX is an attenuating factor for ExB probe data, it is a re-distributing factor for 
Faraday probe data. This is because CEX ions have a different trajectory than the ion it is born from and 
will tend to show up in the wings of a Faraday probe trace. The region considered the “wings” of a 
Faraday probe trace can be seen in Fig. 7 as the region outside of the two outermost dashed vertical 
lines. An ExB probe does not see CEX ions because of a combination of small acceptance angle and a 
tendency to de-emphasize low speed ions. The ExB filter resolution has a dependence on the square of 
the particle speed and CEX ions have only a small fraction of the speed of beam ions. The best way to 
correct for the re-distributing behavior of CEX on Faraday probe traces is to take data at multiple 
distances and facility background pressures. 25 This could not be carried out for the current study due to 
time constraints. Instead, an extrapolated-wing approach is used to remove the parts of the trace that are 
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believed to be dominated by re-distributed CEX ions. 
The extrapolated-wing approach begins by plotting the 
ion current density on a semi-log scale. Curve-fits are 
then performed on the data from a polar angle of 10° 
to 25°. The curve-fit result is then extrapolated to 90°. 
Figure 7 shows a graphical representation of the 
extrapolated-wing method. The black solid line is the 
raw data. The dashed vertical line in the middle is the 
center of the trace (0 is not quite equal to 0° here due 
to a minor mis-alignment in the physical setup). The 
left two and right two dashed vertical lines indicate the 
boundary for the data used to perform the curve-fit on 
the respective side. The light green solid lines 
represent the curve-fit results. 



Polar angle, deg. 

Figure 7. Illustration of the extrapolated-wing 
method for analyzing Faraday probe data. 


In past studies, low secondary electron 
emission (SEE) yield materials like molybdenum and 

tungsten have been used for the collector. 7 ’ 26 ’ 27 Table 1. Summary of SEE data for xenon 

Secondary electrons born on a negatively biased probe io n bombardment of molybdenum. 

will accelerate away from the probe. This effect adds extra 
current to the probe measurement that is indistinguishable 
from the collected ion current. By selecting a low SEE 
yield material, researchers had believed that the effect of 
SEE is negligible and can be ignored. While singly- 
charged xenon-induced SEE yield for molybdenum and 
tungsten are negligible (0.013 to 0.022), the doubly-charged xenon-induced SEE yield are roughly 10 
times higher than the singly-charged yield, and the triply-charged SEE yield are roughly 35 times higher 
than the singly-charged yield. 28 ' 0 Furthermore, Hagstrum discovered that metastable singly-charged 
xenon induces roughly the same SEE yield as doubly-charged xenon. 31 Thus, as long as the plume of the 
thruster is composed of mostly ground-state singly-charged xenon ions, the assumption of negligible SEE 
effect holds. 


Bombarding SEE yield of 


particle 

molybdenum 

Xe + 

0.022 

Xe 2+ 

0.20 

Xe 3+ 

0.70 


During the present testing, the assumption of negligible SEE effect did not hold. The material 
used for the present study is molybdenum, which is considered a low SEE yield material. The reason 
SEE effects are not negligible is because high-power Hall thrusters appear to produce 20-30% doubly- 
plus triply-charged ions by current fraction. Combined with the high SEE yield associated with 
bombardment by higher charge ions, the effect on Faraday probe measurements is in the range of 5- 
10%. To correct for SEE effects on the Faraday probe measurement, we turn to data published by 
Hagstrum. Table 1 summarizes the SEE yield values used in the data analysis of the present study. The 
singly-charged and doubly-charged xenon-induced yields are averages of the SEE yield data for ion 
energies in the range of 200 to 800 eV in Hagstrum’s 1956 work on molybdenum. 29 For both of these 
parameters, the value measured by Hagstrum varied by no more than 10% of the listed average. A 
published value for the triply-charged xenon-induced yield of molybdenum could not be found. The value 
in Table 1 is a projected value based on the similarity in yield between tungsten and molybdenum. The 
ratio of triply-charged induced yield to doubly-charged induced yield for tungsten is 3.5, so the yield for 
molybdenum is projected to be 3.5 * 0.2, or O.7. 28 The search for a more accurate triply-charged xenon- 
induced SEE yield of molybdenum is left for future work. 

Equation (27) shows the relationship between the actual ion current density, the various SEE 
yield values, and the ion current density measured by the Faraday probe. 


J FP 

1 + Yeff 



(27) 


where J is the actual ion current density, J FP is the current density measured by a nude Faraday probe, 
Yeff is the effective SEE yield, and y k is the SEE yield associated with bombardment by the k-th species. 
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Using Eq. (27), we can write a new equation to show the definition of the effective SEE yield for the 
Faraday probe, shown in Eq. (28). 


V 1+ Yk 

For the present study, the effective SEE yield range from 0.06 to 0.09, which is not a negligible amount. 
The amount of metastable singly-charged ions is assumed to be negligible. This assumption is an 
extrapolation from controlled ion beam work done by Hagstrum. 31 In his experiment, Hagstrum varied the 
energy of the electrons used to generate his singly-charged ion beam from 10 to 70 eV. The SEE yield he 
measured jumped quickly up from 0.022 to 0.025 when the electron energy was ramped up from 25 to 30 
eV, and then plateaued out at 0.025 all the way up to 70 eV. While the Hall thruster ionization zone and 
an ion beam discharge chamber are not exactly the same they do share many of the same operating 
principles, including a reliance on impact bombardment ionization. Since the amount of metastable ions 
was so small that the SEE yield increased by only 0.003 for Hagstrum’s experiment, we assumed the 
effect will be of similar magnitude in a Hall thruster and is therefore negligible. 

F. Phenomenological Efficiency Model 


A phenomenological efficiency model is a way of breaking down and describing various physical 
phenomena that affect the overall efficiency of a thruster. The purpose of a phenomenological efficiency 
model is to help researchers identify trends in the processes that drive the overall thruster efficiency. 
Many efficiency models have been proposed for the Hall thruster in the past. The complexity of these 
models depended on the operating environment and the state of knowledge in the community at the time. 
The model used in this paper is the same as a prior work by Shastry, 9 which has evolved over time from a 
number other studies. 8, 4,24,32 It is shown in Eq. (29) to (34). 
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(30) 

(31) 

(32) 

(33) 

(34) 


where r\ a is the anode efficiency, which can also be calculated via Eq. (1), q v is the voltage utilization 
efficiency, r) d is the divergence efficiency, r| b is current utilization efficiency, q m is the mass utilization 
efficiency, r| q is the charge utilization efficiency, V RPA is the ion energy per charge measured by the RPA, 
l b is the total beam current measured by the Faraday probe, m Xe is the mass of xenon, and a m is the part 
of mass utilization efficiency that is calculated from ExB probe data. The voltage utilization efficiency 
describes, on the average, how much of the voltage drop provided by the discharge supply is actually 
used to accelerate the ions. The divergence efficiency describes how much of the kinetic energy imparted 
to the ions is axial, thrust-producing, kinetic energy. The current utilization efficiency describes how much 
of the discharge current is carried by ions instead of electrons. Electrons generate negligible thrust 
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compared to the ions. The mass utilization efficiency describes how much of the mass flow exiting the 
thruster channel is in the form of ions. The charge utilization efficiency is a number of terms representing 
the effects of having multiply-charged species that are not already described by the other terms. 

G. Uncertainty Analysis 

Due to time constraints, a full uncertainty analysis will not be carried out for every element of this 
study. Only uncertainties associated with novel parts of the study and uncertainties that are dominant will 
be analyzed. The uncertainty in the anode efficiency as calculated by Eq. (1) from the measured thrust is 
dominated by the uncertainty in the thrust and is ~4% when calculated via standard error propagation. 

The uncertainty in the voltage utilization as calculated from Eq. (30) is ~1%. This uncertainty is a 
combination of the uncertainty in the RPA analysis method and the Langmuir probe analysis method used 
to correct the RPA measurement. The uncertainty in the charge utilization as calculated from Eq. (34) is 
expected to be less than 0.5%. That is because this value is very insensitive to the value of the current 
fraction. For the current study, the charge utilization value varied from 0.974 to 0.98. 

The uncertainty in the divergence efficiency, current utilization efficiency, and mass utilization 
efficiency varies greatly from condition to condition. Specifically, the ExB probe trace for the 200 V, 50 A, 
300 V 50 A, and 300 V, 67 A cases are highly abnormal. The abnormality will be shown and discussed in 
the results section. With the exception of these three cases, the uncertainties of the aforementioned three 
efficiencies are as follows. 

For the divergence efficiency, the important sources of uncertainty include background CEX 
effects and SEE effects. Based on prior work by Brown, Faraday probe data taken at 5 MTD away from 
the thruster should have an uncertainty of no more than 2% for the facility background pressure found 
during the current study. 25 However, Brown’s work was done on a lower power thruster (nominal 6 kW) 
compared to the 300M (20 kW) and with a correspondingly shorter distance between the probe and the 
thruster. The uncertainty associated with CEX is estimated to be ~3% based on the extra distance 
between the probe and thruster. The effect of SEE on the uncertainty of the divergence efficiency is more 
difficult to calculate due to lack of prior experience. However, as long as the current fraction only 
undergoes small changes as a function of polar angle, the overall error will be small. This is because SEE 
represents a relatively small correction to the local current density. Small changes to a small correction 
should be negligible. Thus, the overall uncertainty in the divergence efficiency is estimated to be 3%. 

Before continuing on to the current utilization, the uncertainties associated with the ExB probe will 
be examined. Based on the equations derived in the appendix, the average velocity resolution of the ExB 
probe is 2.7%. This value is higher for ions with higher energy per charge and is no more than 3.5% for 
the detectable ion populations with the highest energy per charge. This resolution is sufficiently fine so 
that the previously mentioned ExB probe analysis methods are applicable. The ExB probe resolution 
primarily affects broadening of the measured VDF and has minimal effect on the current fractions. The 
uncertainty in the current fraction is expected to be dominated by uncertainty in the average background 
neutral density measured by the ion gauges. The ion gauges are nominally 6% accurate. However, they 
only measure the neutral density in the immediate vicinity. Factors like conductance losses and un- 
compensated temperature effects on the electronics are expected to raise the uncertainty of the neutral 
density measurement to ~15%. Substituting in a ±15% neutral density uncertainty into the CEX correction 
formula yields average uncertainties of ±0.035, ±0.023, and ±0.013 for the Xe + , Xe 2+ , and Xe 3+ current 
fractions, respectively. 

The uncertainty in the current utilization is primarily influenced by uncertainties associated with 
calculation of the effective collector area, CEX effects, and SEE effects. The effective collector area is 
expected to be accurate to ~1% based on different ways in which the gap current can be distributed 
between the collector and the guard ring. The use of the extrapolated-wing approach should remove 
much of the uncertainty associated with the CEX effect, but without additional knowledge about how 
much uncertainty is removed, the same 3% that was used in the calculation of the uncertainty of the 
divergence efficiency will be used here. The uncertainty in the SEE correction depends on the uncertainty 
of the current fractions calculated from the ExB probe data. Substituing in uncertainties in the current 
fractions calculated above into the SEE correction equation (Eq. (27)) the uncertainty in the SEE yield is 
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±0.014. Thus the uncertainty in the current utilization due to SEE correction is ~1.4%. The total 
uncertainty in the current utilization, calculated from standard error propagation, is ~3.5%. 

The uncertainty in the mass utilization is dominated by the uncertainty in the measured anode 
mass flow rate, the total beam current, and the current fraction-dependent factor, a m . The mass flow rate 
measurement has an instrumentation uncertainty of ~1%. The uncertainty in the total beam current is the 
same as the uncertainty in the current utilization (-3.5%). The uncertainty in a m , which is calculated from 
the current fractions, is dominated by uncertainty in the average background pressure. Using ±15% 
neutral density and carrying through the calculations for a m yield an average uncertainty of -2%. The total 
uncertainty, calculated from standard error propagation, in the mass utilization is ~4%. 

Using standard error propagation, the total uncertainty in r\ a calculated via Eq. (29) for operating 
conditions other than 200 V, 50 A, 300 V 50 A, and 300 V, 67 A is ~6%. As previously mentioned the ExB 
probe data for those three operating conditions were abnormal and will be discussed separately. All 
calculations made in this section assume the sources of error are random in nature. 

IV. RESULTS 


A. ExB probe results and comparison between various ExB probe analysis methods 


For the ExB probe analysis in this paper, the new integration approach with the twin Gaussian 
form is assumed to be the most accurate approach, and the associated results will be referred to as the 
baseline results. Results from applying both the new and old (traditional) integration approaches with 
various VDF forms, as well as with the simple peak approach, will be compared to the new integration 
approach with twin Gaussian form. For all ExB probe analysis methods, the analysis program begins by 
finding the location of the four peaks representing the first four charged states. For each charged state 
the left integration boundary is set to one third the distance from the peak being analyzed to the adjacent 
peak; the right integration boundary is set to 1.5 times the peak location. Several sample plots of how the 
ExB probe analysis program functions will be examined next. 

Figure 8 shows a set of analysis plots generated by the program that performs ExB probe 
analysis. The analysis method used is the new integration approach with twin Gaussian form. The 
operating condition is 500 V, 40 A. For convenience, the paper will refer to the six subplots in order from 
left to right, top to bottom as (a) to (f). Subplot (a) shows the raw ExB probe data as black data points with 
red dashed vertical lines showing the location of the first four peaks. Subplots (b), (c), (e), and (f) show 
the four steps in the curve-fitting process. The curve-fitting process fits the appropriate form equation to 
the 4 th peak, subtracts the fitted data from the raw trace, then fits the 3 rd peak, and so on. The data prior 




(b) \ I 

K 

V 

0 50 100 

ExB bias, V 

x 10 -, ° 

(e) fj l 

i\k 

1 

1 

1 

1 

1 


0 50 100 

ExB bias, V 




Figure 8. ExB probe analysis plots for the new integration approach with twin Gaussian form. 

The 300M was operating at 500 V, 40 A. 
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to the fit at each curve-fit step is shown in black dots, the red solid line shows the curve-fit, the magenta 
dashed vertical lines show the curve-fit boundaries, and the blue dashed line shows the leftover result 
after subtraction. Subplot (d) shows the raw data in black dots with a red solid line that represents the 
sum of all four curve-fit steps superimposed on top. Figures 9 to 1 1 show similar plots for the ExB probe 
analysis using the new integration approach with the triangle form, Gaussian form, and variable exponent 
form, respectively. The operating condition is 500 V, 40 A for these three figures. 
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Figure 9. ExB probe analysis plots for the new integration approach with triangle form. 

The 300M was operating at 500 V, 40 A. 








Figure 10. ExB probe analysis plots for the new integration approach with Gaussian form. 

The 300M was operating at 500 V, 40 A. 


16 








Figure 11. ExB probe analysis plots for the new integration approach with variable exponent form. 

The 300M was operating at 500 V, 40 A. 


By comparing subplot (d) of Figs. 8 to 1 1 , the twin Gaussian form appears to provide the best fit 
to the ExB probe data. The variable exponent form also provides a very good fit to the data. Note that a 
good fit to the data does not necessarily imply accurate results, though a bad fit to the data makes the 
analysis results questionable. Having established that the ExB probe analysis program is performing the 
curve-fits correctly, the results from using various analysis methods will now be compared. The parameter 
used for this comparison is a m . This parameter is known to vary the greatest of all the terms in the 
phenomenological model used in the current study. 9 For comparison, data from conditions other than 200 
V, 50 A, 300 V, 50 A, and 300 V, 67 A will be used. Those three conditions will be treated separately 
later. The parameter a m is calculated for each method and for each of the seven operating conditions. 
The calculated values are then subtracted by the baseline value and the differences are divided by the 
baseline value to obtain relative percent differences. Both the average and the standard deviation of this 
relative difference are listed in Table 2. If the average is close to 0, the method in question yields results 
that are generally the same as the baseline. If the standard deviation is close to 0, the method in question 
yields results that closely match the case-by-case variation in the baseline. Thus, for a method to truly 
match the baseline, both the average and the standard deviation of the difference should be close to 0. 

Table 2. Comparison of the relative difference in the a m parameter for 

various ExB probe analysis methods. 

Standard 


New Variable exponent Yes 1 .27% 0.57% 

Old Twin Gaussian Yes 1.80% 0.50% 

Old Triangle Yes 0.38% 0.71% 

Old Gaussian Yes 0.52% 0.72% 

Old Variable exponent Yes 1.21% 0.77% 

Peak height Yes 5.66% 0.92% 

New Twin Gaussian No -17.0% 4.5% 

- Peak height No -7.1% 4.7% 


this division has the effect of cutting down the variation in a m due to variation in the current fraction of 
higher charge states. The low sensitivity of the a m parameter to current fraction is likely a contributing 


Average deviation 


Integration 

approach 

VDF form 

Corrected 
for CEX 

relative 

difference 

of relative 
difference 

New 

Twin Gaussian 

Yes 

Baseline 

- 

New 

Triangle 

Yes 

-0.26% 

0.64% 


New Gaussian Yes -0.51% 0.60% 


A number of 
interesting observations 
can be made from the 
results shown in Table 2. 
First, the difference 
between using the triangle, 
Gaussian, variable 

exponent, and the twin 
Gaussian forms is very 
small. Given the 
uncertainty in the probe 
data, there is no 
meaningful difference in a m 
for the different VDF forms 
used. Recall the a m 
parameter involves dividing 
current associated with 
higher charged states by 
charge. Mathematically, 
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factor to the lack of meaningful differences in the choice of VDF forms. Another possible explanation for 
the small variation is that the exact equation form is not very important as long as the height and width of 
the distribution is properly captured. This explanation is supported by the fact that the triangle form gives 
the closest match to the baseline as the triangle form is only supposed to match the rough height and 
width of the particle distribution. 

Second, the difference between the new and old integration approach is also small though it is 
systematic in nature. The old integration approach gives roughly 1% higher a m than the new integration 
approach. Assuming the new integration approach is the more accurate approach, the old integration 
approach overestimates the charged species term in the mass utilization by ~1%. The standard 
deviations in the error for all methods are also fairly small, which indicate that the various methods follow 
the same trend as the baseline. The only method that deviates by a meaningful amount from the baseline 
is the peak height approach, which overestimates a m by 5-6%. This supports the notion that the most 
important aspect of ExB probe analysis is to properly capture the height and width of each peak. 

Third, there is a large difference when CEX is not accounted for. This result highlights the 
importance of CEX correction for high-power Hall thruster tests where the amount of neutrals between the 
probe and the thruster can be quite large. 

The small difference between the new and old integration approach may seem surprising at first 
given the differences in the form of the integration equations in Eq. (20) and (25). The reason for the 
similarity lies in the fundamental scaling of these two integral equations. A quick analysis of this scaling 
using the ion current fraction integral will now be performed. The current fraction calculated from the old 
(traditional) approach will be labeled with 0° and the current fraction calculated from the new approach 
will be labeled with Qj N . Assuming the incoming ion beam is composed of species with very narrow VDFs 
and all charged species undergo the same voltage drop, Eq. (20) can be simplified as shown in Eq. (35), 
and Eq. (25) can be simplified as shown in Eq. (36). 
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The simplified expressions for the two integration approaches are equal. The voltage drop is typically very 
similar between the charged species, thus for plumes with very narrow VDFs, the two approaches yield 
the same answer. Additionally, Eqs. (35) and (36) imply that any difference in the results between the 
new and old approaches arise mostly from differences in the way the two approaches account for the 
spread in the VDF. 

As previously mentioned, the ExB probe traces for the three low-voltage, high-current conditions 
appear abnormal. Figures 12 to 14 show the raw ExB probe traces for the 200 V, 50 A, 300 V, 50 A, and 
300 V, 67 A operating conditions, respectively. Also shown in the figures are red dashed vertical lines that 
indicate the approximate location where various charged species peaks should have appeared in the 
trace. These three traces exhibit strange characteristics where peaks do not always appear where they 
should have and higher-charged state peaks tend to be eclipsed by their neighbors. Analysis was carried 
out on these traces using the new integration approach with twin Gaussian form. However, the 
uncertainty is expected to be particularly large. 
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Table 3 summarizes the ExB probe analysis results. These results use the new integration 
approach with twin Gaussian form and are corrected for CEX effects. The results shown in the table 
confirm the original assessment that for the current data set, the values of the parameters derived from 
ExB probe data are largely insensitive to variations in the current fractions. However, this observation 
comes with the caveat that the ExB probe traces for low voltage, high current conditions appear abnormal 
and have high uncertainties. Conditions with abnormal ExB traces are highlighted with asterisks under the 
discharge voltage column. 



ExB bias, V 

Figure 12. Raw ExB probe trace for the 200 V, 
50 A operating condition. 



ExB bias, V 

Figure 13. Raw ExB probe trace for the 300 V, 
50 A operating condition. 



ExB bias, V 

Figure 14. Raw ExB probe trace for the 300 V, 
67 A operating condition. 


Table 3. Summary of ExB probe data analysis results. 

Asterisks in the discharge voltage column indicate conditions with abnormal ExB traces. 


Discharge 
voltage, V 

Discharge 
current, A 

A, 

Xe + 

Xe 2+ 

O, 

Xe 3+ 

Species 

fraction, 

Xe + 

Species 

fraction, 

Xe 2+ 

Species 

fraction, 

Xe 3+ 

1c 

dm 

*200.6 

49.8 

0.758 

0.224 

0.019 

0.893 

0.103 

0.005 

0.980 

0.876 

300.2 

33.5 

0.812 

0.130 

0.058 

0.931 

0.055 

0.014 

0.980 

0.896 

*300.4 

50.5 

0.714 

0.220 

0.066 

0.879 

0.104 

0.017 

0.974 

0.846 

*301.3 

67.8 

0.738 

0.200 

0.062 

0.890 

0.094 

0.016 

0.975 

0.859 

400.1 

25.4 

0.788 

0.140 

0.072 

0.921 

0.062 

0.018 

0.978 

0.882 

400.7 

37.5 

0.776 

0.139 

0.085 

0.916 

0.062 

0.021 

0.976 

0.874 

400.6 

49.8 

0.785 

0.140 

0.075 

0.919 

0.062 

0.019 

0.977 

0.880 

500.3 

20.2 

0.752 

0.177 

0.071 

0.901 

0.082 

0.018 

0.976 

0.864 

500.7 

30.2 

0.793 

0.145 

0.062 

0.921 

0.064 

0.015 

0.979 

0.886 

500.4 

40.3 

0.759 

0.163 

0.078 

0.906 

0.074 

0.020 

0.975 

0.867 
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B. Faraday probe results 


Review of the Faraday probe results will start with data plots from a few cases. Figures 15 to 18 
show these results. As before, the dashed vertical line in the middle of the plot is the center of the trace. 
The left two and right two dashed vertical lines indicate the data used to perform the curve-fit on the 
respective side. The light green solid lines represent the curve-fit results. 



Polar angle, deg. 

Figure 15. Faraday probe analysis plot, 200 
V, 50 A. 



Polar angle, deg. 

Figure 16. Faraday probe analysis plot, 300 
V, 67 A. 




Polar angle, deg. 

Figure 17. Faraday probe analysis plot, 400 
V, 50 A. 


Polar angle, deg. 

Figure 18. Faraday probe analysis plot, 500 
V, 40 A. 


Generally, the central portion of the Faraday probe trace tends to be more of a sharp peak for 
higher-discharge-voltage cases. This trend indicates the plume is more collimated when the thruster is 
operating at higher discharge voltages. 

The CEX wings appear fairly broad for all cases. Brown previously performed a study of Faraday 
probe traces as a function of background pressure. In that study, the CEX wings can be seen to become 
more prominent with rising background pressure. 25 Although the current study was performed with a 
background pressure that is typically considered low in Hall thruster testing, the Faraday probe data 
shows that the extra distance required for interrogating higher power Hall thruster gave rise to prominent 
CEX wings. Recall that the minimum distance between the Faraday probe and the thruster is limited by 
the need to approximate the source of the plume as a point source. In future high-power Hall thruster 
tests, Faraday probe sweeps at multiple distances may be necessary to obtain accurate data. 
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Table 4. Summary of Faraday probe data analysis results. 

Asterisks in the discharge voltage colmin indicate conditions with abnonnal ExB traces. 


Disch. 

voltage, 

V 

Disch. 

current, 

A 

Original 
Div. Beam 

angle, current, 
deg. A 

Ex.-wing 
Div. Beam 

angle, current, 
deg. A 

Effective 

SEE 

yield 

Beam 
current, 
corrected 
for SEE, A 

Hd 

lb 

* 200.6 

49.8 

30.7 

47.1 

26.3 

42.7 

0.065 

40.1 

0.804 

0.806 

300.2 

33.5 

25.6 

32.2 

15.1 

27.4 

0.067 

25.7 

0.932 

0.766 

* 300.4 

50.5 

27.0 

48.9 

18.4 

41.6 

0.086 

38.3 

0.901 

0.759 

* 301.3 

67.8 

27.9 

66.1 

20.4 

56.8 

0.081 

52.6 

0.878 

0.776 

400.1 

25.4 

25.7 

24.9 

14.9 

21.2 

0.075 

19.7 

0.934 

0.777 

400.7 

37.5 

25.1 

36.7 

14.6 

31.5 

0.081 

29.1 

0.937 

0.776 

400.6 

49.8 

24.6 

48.1 

14.4 

41.2 

0.077 

38.2 

0.938 

0.768 

500.3 

20.2 

25.1 

20.1 

15.3 

17.4 

0.081 

16.1 

0.931 

0.795 

500.7 

30.2 

25.7 

29.7 

15.2 

25.6 

0.072 

23.9 

0.931 

0.792 

500.4 

40.3 

25.1 

39.6 

14.7 

34.0 

0.082 

31.4 

0.936 

0.780 


Table 4 summarizes the Faraday probe analysis results. Divergence angle and beam current 
results are shown for applying versus not applying the extrapolated-wing method. All results beyond the 
6 th column are derived from applying the extrapolated-wing method. As before, operating conditions with 
abnormal ExB traces are highlighted with asterisks under the discharge voltage column. The use of the 
extrapolated-wing approach appears to greatly change both the divergence angle and the beam current. 
The implication is that the CEX wings are broad as previously seen in the Faraday probe plots. Effective 
SEE yield varies from 0.06 to 0.09 (possibly more since ExB probe data for some conditions are 
questionable). 

C. Efficiency analysis results 


The full physical picture of the NASA-300M will now be studied. Table 5 shows the thrust and the 
various efficiencies that were calculated during the current study. The probe-derived n a refers to the 
anode efficiency derived from the use of probe data with the phenomenological model (Eqs. (29) to (34)). 
The thrust-derived r) a refers to the anode efficiency calculated from thrust measurement via Eq. (1). 

Table 5 shows several signs that indicate the abnormality in the ExB probe traces for the low- 
voltage, high-current operating conditions is problematic. The mass utilization efficiencies for the 200 V, 
50 A and 300 V, 67 A are well above 1 and unphysical; one of the possible causes of these results is the 
inability of the ExB probe analysis program to properly isolate the peaks associated with the higher 
charged states, which causes an overestimation of the a m parameter and the q m . Another possible culprit 
is the broad CEX wings present in the Faraday probe traces for these operating conditions. Comparing 
Figs. 15 and 16 to Figs. 17 and 18, the CEX wings have much higher amplitude relative to the central 


Table 5. Summary of thrust measurement and efficiency analysis. 

Asterisks in the discharge voltage column indicate conditions with abnormal ExB traces. 


Disch. 

voltage, 

V 

Disch. 

current, 

A 

Thrust, 

mN 

Hv 

Hd 

Hb 

Hm 

He 

Probe- 

derived 

Ha 

Thrust- 

derived 

Ha 

Difference 

* 200.6 

49.8 

671 

0.873 

0.804 

0.806 

1.198 

0.980 

0.665 

0.565 

0.100 

300.2 

33.5 

652 

0.926 

0.932 

0.766 

0.969 

0.980 

0.628 

0.654 

- 0.026 

* 300.4 

50.5 

912 

0.934 

0.901 

0.759 

1.050 

0.974 

0.652 

0.651 

0.001 

* 301.3 

67.8 

1161 

0.906 

0.878 

0.776 

1.178 

0.975 

0.709 

0.632 

0.077 

400.1 

25.4 

581 

0.942 

0.934 

0.777 

0.938 

0.978 

0.627 

0.657 

- 0.031 

400.7 

37.5 

852 

0.940 

0.937 

0.776 

0.991 

0.976 

0.660 

0.692 

- 0.031 

400.6 

49.8 

1099 

0.943 

0.938 

0.768 

1.049 

0.977 

0.696 

0.693 

0.003 

500.3 

20.2 

528 

0.962 

0.931 

0.795 

0.913 

0.976 

0.634 

0.668 

- 0.034 

500.7 

30.2 

770 

0.942 

0.931 

0.792 

0.990 

0.979 

0.673 

0.673 

0.000 

500.4 

40.3 

1004 

0.959 

0.936 

0.780 

1.018 

0.975 

0.695 

0.686 

0.009 
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peak height for the operating conditions with abnormal ExB probe traces than for those with normal 
traces. The analyzed data may still contain a non-negligible amount of CEX broadening despite the use of 
the extrapolated-wing method. The mass utilization efficiency calculated for certain other operating 
conditions also rises above 1 though they are still physical if one accounts for the uncertainty in the 
measurement. The fact that the mass utilization efficiency is generally very high may be an indication that 
the CEX effect on the Faraday probe data is not completely removed. No attempt is made to artificially 
reduce the mass utilization efficiencies that exceeded 1. If the error sources are truly random, as was 
assumed in the uncertainty analysis, values that rise above 1 due to random error in one efficiency may 
be offset by values that are below average in another efficiency. Artificial pruning of the efficiency values 
can create extra systematic error. As a result of the abnormally high mass utilization efficiency for the 200 
V, 50 A and 300 V, 67 A cases, the probe-derived anode efficiency is also abnormally high for those two 
cases. 


Setting aside the low-voltage, high-current operating conditions, we can study Table 5 for trends. 
Recall the original purpose of the phenomenological model is to provide a means of studying trends and 
not necessarily to match the thrust-derived anode efficiency. Table 5 shows several trends that have 
been previously observed for low and medium power Hall thrusters. First, the voltage utilization tends to 
be higher for higher discharge voltage operation. This is because the voltage needed to extract electrons 
from the cathode does not vary greatly with the discharge voltage. This extraction voltage takes up an 
ever smaller fraction of the discharge voltage with increasing discharge voltage. 

Second, plume divergence decreases with increasing discharge voltage. Whatever radial voltage 
drop that is experienced by the ions is likely not varying much with discharge voltage. As a result, the 
axial voltage drop become more dominant as the discharge voltage increases and the beam becomes 
more collimated. 

Third, mass utilization efficiency appears to rise with mass flow rate. This is expected because 
the density of the Hall current increases with the discharge current, which increases with the mass flow 
rate. With that said, given the previously described issues with parts of the ExB probe data and CEX 
effects on the Faraday probe data, some part of the observed trend may be an artifact of one or both 
issues. 


Fourth, current utilization efficiency and charge utilization efficiency are both fairly constant 
across all operating conditions. Although this trend is not absolutely certain due to CEX-induced 
uncertainties, the insensitivity of the charge utilization efficiency to operating condition has been 
previously observed. 9 

Although the main purpose of the phenomenological model is not to reproduce the thrust-derived 
anode efficiency, comparison of the anode efficiency provides a good reality check on the results. The 
last column in Table 5 shows the difference between the probe-derived and the thrust-derived anode 
efficiencies. With the exception of the operating conditions with abnormal ExB data, the two anode 
efficiencies agree to within ~4%. Recall the total uncertainty in the thrust-derived anode efficiency is ~4% 
and the total uncertainty in the probe-derived anode efficiency is ~6%. 

V. CONCLUSION 

A far-field plume characterization study has been carried out on the NASA-300M. In-depth work 
on the analysis methods was carried out to identify factors that may have been previously ignored but are 
important when analyzing the plume of a high-power Hall thruster. 

Fresh derivation of the formulas that describe the operation of an ExB probe in the plume of a 
Hall thruster revealed a new integration approach for calculating current fraction and density fraction. 
Further study showed that the new integration approach yields the same answer as the old integration 
approach for a mono-energetic beam. Differences in results between the two approaches are attributed to 
the different ways the two approaches treat the spread in the particle VDF. Choice of different VDF forms 
had little effect on the overall results, suggesting that the result is accurate as long as the height and 
width of the particle distributions are properly captured. While the two integration approaches agree fairly 


22 



well with each other, they disagree with the peak height approach. Additional analysis show that the CEX 
correction factor made the biggest difference in the results. 

The plume of the NASA-300M is found to carry fairly sizable doubly- and triply-charged 
populations. Since these higher-charge-state ions have non-negliaible SEE yield, a new SEE correction 
formula for nude Faraday probe is derived from data by Hagstrum. 31 This formula predicted an effective 
SEE yield of 0.06 to 0.09 for the conditions present in the current study. 

The trends revealed by the phenomenological model generally agree with past observation on 
lower power thrusters. When the thruster operates at higher discharge voltages, it tends to utilize more of 
the discharge voltage and have lower divergence. When the thruster operates at higher mass flow rate, it 
tends to utilize more of the mass flow. Barring a few conditions with abnormal data, the anode efficiency 
derived from probe measurements agree with that derived from thrust measurements to within 4%. 

From the present study, one can conclude that for high-power Flail thruster plume 
characterization, extra care must be taken with positioning of the probe. Unlike low-power Flail thruster 
studies, thruster heating places a demanding constraint on how close to the thruster the ExB probe and 
RPA can be placed. Flowever, if the probes are too far away, CEX effects can degrade the quality of the 
probe data (this may have been the cause of the abnormal ExB probe traces). Probe shields, shutters, 
and possibly even active cooling may become necessary for characterization studies of very-high-power 
Hall thrusters (50+ kW). For the Faraday probe, the need for the point-source approximation will drive the 
closest distance to the thruster that the probe can be used. Given the large influence of CEX wings on the 
analysis results, it is likely that scanning at multiple distances and background pressures may become 
necessary for obtaining accurate results. 


VI. FUTURE WORK 

Further work is needed in several areas. Now that the AFRL RPA has been validated, the far-field 
tower is no longer limited by the usable range of the GRC RPA. Future testing can involve putting the 
tower closer to the thruster. The ExB probe data may be of higher quality when taken with the probe 
closer to the thruster. 

Given the results of the current study, future work should involve taking Faraday probe data at 
multiple distances, possibly with multiple background pressures. However, if CEX effects continue to 
overwhelm the wings of the Faraday probe traces, a gridded Faraday probe may be needed in 
conjunction with a nude Faraday probe. Such a setup has the added benefit of being able to suppress 
SEE. For SEE correction of nude Faraday probe data, more research must be carried out to determine a 
proper value of triply-charged xenon-induced SEE yield. 

The next far-field plume characterization test will be carried out on the NASA-300MS, the 
magnetically shielded version of the NASA-300M. While some of the above changes will be implemented 
for this next test, others may not be implemented in order to make the results between the two tests 
comparable. 
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APPENDIX 


Derivation of ExB probe resolution begins with 
an assumption on the geometric design of the ExB 
probe. Figure 19 shows the design of a generic ExB 
probe with a drift section. This simplified and 
exaggerated diagram also represents the ExB probe 
used for the current study. On the left of the diagram, 
there is the collimator, which has a circular entrance 
orifice of diameter a u a circular exit orifice of diameter 
a 2 , and a length of L c . The particle is injected at the 
maximum collimator acceptance angle, a, and enters 
the filter section. The filter section is of length L f and is 
sufficiently wide that the particle does not hit the bias 
plates (grey rectangles). The electric and magnetic 
fields in the particle path are uniform, perpendicular, 
and of magnitude E 0 and B 0 , respectively. The example particle in this derivation is traveling slightly faster 
than the filter velocity u 0 , which is equal to E 0 /B 0 . The velocity of the particle is assigned to be u 0 +Au. 
Because the particle is slightly too fast, it experiences a net downward acceleration in the filter section 
and passes into the drift section. The drift section is of length L d . The particle is just slow enough to pass 
through the circular orifice of diameter a 3 in the SEE suppression plate. The collector is sufficiently large 
as to collect all ions passing through the orifice on the SEE suppression plate. In some ExB probe 
designs, an exit collimator replaces the SEE suppression plate. The diameter of the entrance orifice for 
the exit collimator would become a 3 . Two assumptions will be used in this derivation. One, the 
acceptance angle is small enough that small angle approximation applies. Two, Au « u 0 , so the equation 
can be linearized. 

Applying the equations of motion to the collimator section, we can obtain Eq. (37) for the 
acceptance angle of the collimator. Note the use of the small angle approximation. 
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Figure 19. Simplified diagram of a generic 
ExB probe with drift section. 
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(37) 


Applying the equations of motion to the filter section, we obtain Eq. (38) for the x-axis and Eq. (39) for the 
y-axis. Note the use of the Au « u 0 assumption in Eq. (38). 


Tf ~ Lf /u 0 

Ay f = A Aax f 2 +u 0 sinax f 


Au = — - 


m 

ZeB n 


-Aa 


(38) 

(39) 

(40) 


where x f is the dwell time of the particle in the filter section, Ay f is the change in y-position of the particle 

between where it enters and where it exits the filter section, and Aa is the net acceleration on the particle 

due to the velocity of the particle deviating from the filter velocity. The relationship between Au and Aa is 
derived from force balance in the filter section and shown in Eq. (40). Applying the equations of motion to 
the drift section, we obtain Eq. (41) for the x-axis and Eq. (42) for the y-axis. 

x d «L d /u 0 (41) 

A yd =u y,d T d =Aax f T d +u 0 sinaT d (42) 


where x d is the dwell time of the particle in the drift section, Ay d is the change in y-position of the particle 
between where it enters and where it exits the drift section, and u y , d is the y-velocity of the particle as it 
enters the drift section. u y , d is calculated from the equations of motion for the filter section. Next, we note 
that the total change in the y-position of the particle between where it enters the filter section and where it 
exits the drift section is equal to the negative of the difference in y-position between the upper lip of the 
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collimator exit and the lower lip of the orifice on the SEE suppression plate. This relationship is expressed 
in Eq. (43). 


r 

Ay f +Ay d =- 

V 


a 2 +a 3 i 


2 


J 


(43) 


Next, substitute Eqs. (37)-(39), (41), and (42) into Eq. (43) and re-arrange some terms to obtain Eq. (44). 
Lastly, substituting Eq. (44) into Eq. (40) to obtain an expression for the maximum deviation in velocity 
that a particle can have and still reach the collector. 


Aa _ _ a 2 +a 3 + ( a 1 + a 2)(*-f + Ld)/L c i u ‘ 

L Lf 2 +2L d L f 

mu Q 2 a 2 + a 3 + ( a 1 + a 2)(Lf + l-d)/l-c 

ZeB 0 | L f 2 + 2L d L f 


Au = 


(44) 

(45) 


Note, if the acceptance angle of the ExB probe is larger than the angle subtended by the plasma source, 
the angle a is no longer the limiting parameter. Instead, the angle (3 (Eq. (46)) should be used and the 
form of Eq. (45) changes to Eq. (47). 


sinp « 


D PS 

2 R ExB 


Au = 


mu Q 2 

ZeB 0 


a 2 +a 3 +2sinp(Lf +L d ) 
Lf 2 +2L d L f 


(46) 

(47) 


where D PS is the diameter of the plasma source and R ExB is the distance between the ExB probe and the 
plasma source. For the probe used in this study, the acceptance angle is sufficiently small that Eq. (47) 
does not apply. 

For an ExB probe without a drift section, we can set L d to 0 in Eq. (45), which simplifies into Eq. 

(48). 


Au = 


mu Q 2 

ZeB 0 


a 2 + a 3 + (^1 + a 2) * Lf /L c 


(48) 


By comparing Eq. (48) to Eq. (45), one can see that the addition of the drift section will generally improve 
the resolution of the ExB probe. Note that a u a 2 , a 3 , L c , L f , L d , and R Ex b are all geometric constants for a 
given experimental setup. D PS is also typically a constant (though one can never be too careful when 
generalizing about plasma sources). A scaling equation for the maximum accepted deviation in velocity of 
a generic ExB probe, shown in Eq. (49), can now be written. 


Au = 


mu 0 
Ze Bn 


(49) 


where G is a geometric constant. 

Here are some quick observations about the maximum accepted deviation in velocity for an ExB 
probe. One, the numerator of Eq. (49) contains the ion kinetic energy while the denominator contains the 
ion charge. In other words, Au scales with the ion energy per charge of the incoming particles, which is 
slightly lower than the discharge voltage of a Hall thruster. The implication of this trend for a Hall thruster 
is that Au is species-insensitive because all charged species undergo similar voltage drop and have 
similar ion energy per charge. When calculating ExB probe resolution as a relative percentage, the value 
for the singly-charged species will always be the biggest. That is because Au is roughly constant while u 0 
increases with species charge. 

Two, dividing both sides of Eq. (49) by u 0 , we obtain Eq. (50). 
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( 50 ) 


Au _ mu 0 
u 0 ZeB 0 

This equation says that for a given charged species, the ExB probe resolution in relative percentage 
increases with the velocity of the particle. For instance, an ExB probe with a ±2% resolution when probing 
the plume of a 250-V thruster would have a ±4% resolution when probing the plume of a 1000-V thruster. 
The probe resolution in relative percentage scales as the velocity of the particles, which scale as the 
square root of the accelerating voltage. This is perhaps part of the reason why ExB probe traces for ion 
thruster are never able to resolve the VDF of the plume particles. The high accelerating voltage (1000+ V) 
associated with an ion thruster means Au will generally be large. The combination of low VDF spread and 
large Au means that Au is generally larger than the VDF spread for an ion thruster. This interesting trend 
also means that an ExB probe used to interrogate the plume of a Hall thruster operating at high discharge 
voltage should have tight dimensions to ensure Au is several times smaller than the width of the ion VDF. 
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Introduction (1/2) 



• NASA Space Technology Mission Directorate (STMD) is investing in the 
development of a 15-kW Hall thruster, 

■ Stepping stone toward higher power system for human exploration 

■ Provide propulsion system for Solar Electric Propulsion Technology 
Demonstration Mission (SEP TDM), which demo SotA solar array, electric 
propulsion, and associated spacecraft technologies. 

• NASA GRC and JPL jointly developing high-power magnetically-shielded 
Hall thruster via combined experiment-simulation approach 
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Introduction (2/2) 



• To guide development, a magnetically-shield version of the NASA- 
300M called NASA-300MS will be built and tested 

• Experimental work include extensive far-field, near-field, and 
interior probing. 

• The present work covers the far-field plume characterization of the 
300M, which will serve as a baseline for comparison to the 300MS 
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Outline 



• Introduction 


• Experimental Setup 

• Data Analysis 

■ New integration formula for ExB probe analysis 

■ Effects of SEE for far-field Faraday probe 

• Results 


• Future Works and Conclusion 
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NASA-300M Hall thruster 









NASA-300M is a 20-kW Hall thruster: 

■ Magnetic lens topology, symmetric about 
channel centerline, 

■ Magnetic field strength optimized for anode 
efficiency 

■ Centrally-mounted cathode, 8% CFF 

■ Nominal operating condition: 500 V, 40 A 

10 conditions tested 



■ Discharge power: 10, 15, 20 kW 

■ Discharge voltage: 200 - 500 V 

■ Discharge current: 20 - 67 A 

■ Thrust: 530- 1160 mN 

■ Anode efficiency: 57 - 69 % 
Tested in VF5 

7 ion gauges to monitor pressure 

■ Pressure near thruster < 3e-5 Torr 



Discharge current, A 
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Setup: Facility and motion system 


Farfield probe tower 

,\ V 

r<j\ i 




+R 




.* ** T ■ 

Tri* 



Faraday 

probe 




lY 
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Probes 




' — Insulator 


• Probe tower at ~27 MTD 

■ RPA 1, GRC design, cannot be 
placed closer than 25 MTD 

■ RPA 2, new AFRL design 

■ ExB, commercial 

• FP polarly-swept at 5 MTD 
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Data analysis: 
Langmuir probe and RPA 



• Langmuir probe analysis: simple Langmuir probe theory 

■ 20 traces per test point, averaged and then smoothed 

■ Plotted derivative of probe current with respect to probe voltage 

■ Peak-amplitude voltage is plasma potential 


• RPA analysis: 

■ The trace is smoothed 

■ Plot -dC/dV versus RPA bias 
voltage 

■ Traditionally, take the most 
probable voltage as the ion 
energy per charge 

■ For this study, took ensemble 
average of the part of the 
peak that exceeded 50% of 
max amplitude 

RPA Analysis plot, 200V, 50 A^ 



0 100 200 300 400 500 

RPA bias, V 


x 10 
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Data analysis: 

ExB probe 

• Traditional approach involves fitting the species peaks in the 
collector current versus bias voltage plot to various equation forms 
and calculating the area under the curve 

• New approach is very similar except the integration is carried out 
using new integration formulas that were vigorously derived 

• Steps for analysis of ExB probe data: 

1. Determine rough locations of the peaks 

2. Choose curve-fit boundaries 

► 1/3 of the way to the next peak on left, 1.5X peak location on right 

3. Perform curve-fit using the chosen VDF form 

4. Integrate the fitted curve using the new integration formulas to 
obtain raw current and species fractions 

5. Correct for CEX attenuation 

► CEX model require background neutral density (from ion gauges), 

► and average ion energy per charge (from RPA) 
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Data analysis: 

ExB probe - VDF forms 



Tried four VDF forms 

■ Triangle: 

■ Gaussian: 

■ Variable exponent: 

■ Twin Gaussian: 


fi( u ) = 


(-a/b)*|u-c| + a, c-b<u<c + b 
0, else 


fj(u) = a*exp[-b(u-c) 2 ] 
fj(u) = a*exp[-b|u-c| n ] 

fj (u) = a-| *exp[-b-| (u-c-|) 2 ] + a 2 *exp[-b 2 (u-c 2 ) 2 ] 


Perform curve-fit using this equation (insert fj from above): 

l P>j (Vp) = K 1 n j Vp 3 f,(K 2 V P ) 

Then use these formulas for species and current fractions: 


J »oo 

0 


Ip.i(Vp) 


Ci - 


0 Vp 3 


dVc 


Z: 


poo 

Jo 


Ipj(Vp) 


Z n k y r 

k 4^ Jo 


■p,k(Vp) 

Vp 3 


Qj = 


o V P 2 


dV F 


dV F 


Integrating area under the curve is 
equivalent to: 


Q; = J 


X'k V7 f°°W V P) 

k r Jo v P 2 

/•OO 

[ lp,j(Vp)dVp 


?' k Ej n ”p,k(V P )dVp 


dVc 


*Derivations for all equations shown in paper. 
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ExB current, A ExB current, A 


Data analysis: 

ExB probe - curve fit (1/2) 










Black dot = data point 
Vertical dashed line = 
peak locations 
Red line = sum of 
curve-fits 


Black dot = data point 

Red line = curve fit 

Blue dashed line = left-over 

Pink vertical line = curve-fit boundaries 
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ExB current, A ExB current, A 


Data analysis: 

ExB probe - curve fit (2/2) 



Triangle form 



0 50 100 

ExB bias, V 


Variable exponent form 



Gaussian form 



0 50 100 

ExB bias, V 


Twin Gaussian form 



0 50 100 

ExB bias, V 
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Data analysis: 
Faraday probe 







Assumptions: 

■ Plume is axisymmetric 

■ Constant charge-state ratio 

Charge-weighted divergence angle equation for polar sweep: 


< cos 5 >= 


9 pTt/2 

27iR FP j(0) cos0sin0d0 

9 pTc/2 

27i Rpp J j(0) sin0d0 


• Denominator equal to 
total ion beam current 

• Extrapolated-wing to 
counter effects of CEX 
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Data analysis: 

Faraday probe - SEE effects 



• SEE effects are negligible if plume 
consist mostly of singly-charged 
ions, which is not true for this 
study 

• Data for Xe + , Xe 2+ yield from 
experiments by Hagstrum; Data 
for Xe 3+ yield extrapolated from 
said experiments* 

• Figure from Hagstrum's work* 


J = 


'FP 


= J 


FP 


Si 


Qi 


1 + Yeff " V 1 + Yk 

Equation above used to correct 
SEE effects, 

■ y eff = effective SEE yield for given 
charged species composition; 
0.06 to 0.09 for the present study 


Bombarding 

SEE yield of 

particle 

molybdenum 

Xe + 

0.022 | 

Xe 2+ 

0.20 

Xe 3+ 

0.70 | 



'Hagstrum, Phys. Rev., Vol.104, No. 3, p.672; Hagstrum, Phys. Rev., Vol.96, No. 2, p.325 
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Data analysis: 

Phenomenological efficiency model 

• A phenomenological efficiency model is used to help study trends: 

T 2 

r|a ~ D = Tlq^vTIdTIb^lm 

^ m a K d 



Charge utilization: 



Voltage utilization: 

Plume divergence: 
Current utilization: 

Mass utilization: 


T|v =1 “ 



'Hd - (< COS 8 >) 


2 


= 


Tim = 


]b 

Id 

Mxe’b 

rh a e 



z k 


T = thrust 

m a = anode mass flow rate 
P d = discharge power 

= k-th species current fraction 
Z k = k-th species charge state 
V L = loss voltage 
V d = discharge voltage 
5 = divergence angle 

l b = total ion current 
l d = discharge current 
M xe = mass of xenon atom 
e = elementary charge 
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Results: Abnormal ExB traces 



200 V, 50 A 



300 V, 50 A 



• Higher charge species could 
not be isolated from lower 
charge species 

• Analysis still carried out the 
same way 

• Uncertainties expected to be 
high 


300 V, 67 A 
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Results: ExB results 



• Comparison made using a m , which is a term in r| m -> a m = 

k Zk 

• On the average, new integration approach ~1% less than old 

■ Turns out they are equal for mono-energetic beam (derived in paper) 

■ Difference is due to difference in treatment of VDF spread 

• Different VDF forms gave roughly the same results 

■ Exact shape not important if height and width properly captured 


Integration 

approach 

VDF form 

Corrected 
for CEX 

Average relative 
difference in a m 

Std. dev. of relative 
difference in a m 

New 

Twin Gaussian 

Yes 

Baseline 

- 

New 

Triangle 

Yes 

-0.26% 

0.64% 

New 

Gaussian 

Yes 

-0.51% 

0.60% 

New 

Variable exponent 

Yes 

1 .27% 

0.57% 

Old 

Twin Gaussian 

Yes 

1 .80% 

0.50% 

Old 

Triangle 

Yes 

0.38% 

0.71% 

Old 

Gaussian 

Yes 

0.52% 

0.72% 

Old 

Variable exponent 

Yes 

1.21% 

0.77% 

- 

Peak height 

Yes 

5.66% 

0.92% 

New 

Twin Gaussian 

No 

-17.0% 

4.5% 

- 

Peak height 

No 

- 7.1% 

4.7% 


A Only ExB probe traces that are not abnormal are used 
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Results: Efficiency Analysis 



• Asterisk in discharge voltage column indicate conditions with 
abnormal ExB traces 

■ 2 of 3 shows abnormally high mass utilization (in orange) 

• n v and r| d increase with discharge voltage 

• n m increase with anode mass flow rate 


Disch. 

voltage, 

V 

Disch. 

current, 

A 

Thrust, 

mN 

1v 

n d 

n b 

1m 

n q 

Probe- 

derived 

la 

Thrust- 

derived 

n a 

Difference 

*200.6 

49.8 

671 

0.873 

0.804 

0.806 

1.198 

0.980 

0.665 

0.565 

0.100 

300.2 

33.5 

652 

0.926 

0.932 

0.766 

0.969 

0.980 

0.628 

0.654 

-0.026 

*300.4 

50.5 

912 

0.934 

0.901 

0.759 

1.050 

0.974 

0.652 

0.651 

0.001 

*301.3 

67.8 

1161 

0.906 

0.878 

0.776 

1.178 

0.975 

0.709 

0.632 

0.077 

400.1 

25.4 

581 

0.942 

0.934 

0.777 

0.938 

0.978 

0.627 

0.657 

-0.031 

400.7 

37.5 

852 

0.940 

0.937 

0.776 

0.991 

0.976 

0.660 

0.692 

-0.031 

400.6 

49.8 

1099 

0.943 

0.938 

0.768 

1.049 

0.977 

0.696 

0.693 

0.003 

500.3 

20.2 

528 

0.962 

0.931 

0.795 

0.913 

0.976 

0.634 

0.668 

-0.034 

500.7 

30.2 

770 

0.942 

0.931 

0.792 

0.990 

0.979 

0.673 

0.673 

0.000 

500.4 

40.3 

1004 

0.959 

0.936 

0.780 

1.018 

0.975 

0.695 

0.686 

0.009 
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Future Work 



• Try to take ExB probe data closer to thruster to see if 
the abnormal traces become normal 


• CEX effects on Faraday probe trace may have made 
H m too high, need to take data at multiple distances 
and background pressures 

• Try to find the true value of triply-charged xenon- 
induced SEE yield of molybdenum 

• Perform similar tests on the NASA-300MS 
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Conclusion 



• New derivation of ExB probe analysis revealed a new formula for 
integration 

■ The result is roughly equal to integrating the area under the curve (old 
approach) for the present data set 

• Comparison between four different VDF forms suggest the exact 
form does not matter as long as the height and width are captured 

• Found a formula for offsetting the effects of secondary electron 
emission from a Faraday probe 

■ For this study, 6-9% effect on total ion current 

• Really need to take data at multiple distances and pressure to 
correct for CEX effects on Faraday probe data 

• Overall, this 20-kW Hall thruster have similar physical trends as 
lower power thrusters 

■ n v and n d increase with discharge voltage 

■ n m increase with anode mass flow rate 
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